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At the LHC, many new physics signatures feature the pair-production of massive particles with 
subsequent direct or cascading decays to weakly-interacting particles, such as SUSY scenarios with 
conserved conserved R-parity or _ff — > W{lu)W{iv). We present a set of dimension-less variables 
that can assist the early discovery of processes of this type in conjunction with a set of variables 
with mass dimension that will expedite the characterization of these processes. 



I. INTRODUCTION 

The LHC is currently delivering high energy pp col- 
lisions and in the coming months and years it will be 
exploring the TeV scale. Many of the previously unob- 
served processes that experimenters will be searching for 
involve high transverse momentum Standard Model (SM) 
particles, such as leptons and jets, along with missing 
transverse momentum - a characteristic of a large class 
of models with conserved, discrete quantum numbers, re- 
sulting in new particles being produced in pairs and de- 
caying to an even number of stable, weakly-interacting 
particles that will escape detection. 

In the past years the development of kinematical vari- 
ables that assist the discovery of such processes has been 
intense and rich [T|-|15j. Here we introduce the vari- 
ables Mji and M/j* , whose distributions contain informa- 
tion about the masses of pair-produced particles and the 
weakly-interacting particles resulting from their decays. 
Additionally, we discuss the dimension-less variables R 
and R* that can be used to select events of interest in 
the presence of large backgrounds, in a number of differ- 
ent exclusive and inclusive final states. 



II. Mr 

To define the we consider a simple example. We 
assume two massive particles, Gi and 6*2, are produced 
through a hard partonic subprocess in a hadron-hadron 
collision. Furthermore, we assume that these two parti- 
cles have the same mass, Mq. In the G1G2 rest frame 
(CM frame) the particles Gi and G2 have equal and op- 
posite momentum, with four-vectors which we define as 

p[Gi\ = pi = Mg7cm{1,/3cm} 

p[G2] = P2 = A^g7cm {1, -/3cm}, (1) 

such that {pi + P2Y = s, where s is the usual Mandel- 
stam variable describing the hard partonic subprocess. 
We further assume that each of the two particles Gi de- 
cays as follows: Gi — >■ QiXi- We assume each Qi is a 
stable, mass-less particle that will be visible in our de- 
tector. Each Xi is assumed to be stable, potentially mas- 
sive (with mass M^) and weakly interacting such that it 



escapes detection. In their respective Gi rest frames, the 
decay products of each Gi have four-momenta defined as 

p[Qi\ = qi = 

p[W,] = lo, = ^{Rg^,-u,}, (2) 

where Ma = —37^, Rgx = m^-m^ ^'^'^ ^^^"^ 
a unit vector. To go from the rest frame of Gi (G2) to 
the CM frame, qi and uji ((72 and 0^2) are boosted to a 
frame traveling at a velocity /3cm (—/3cm) with respect 
to the Gi (G2) rest frame. Finally, we assume that to 
move from the CM frame to the lab frame, each of the 
final state particles is boosted to a frame traveling at a 
velocity Pl = {Pt, Pi), where Pt and Pi are the transverse 
and longitudinal components of this boost respectively. 
The transformations taking the final state particles from 
their respective Gi rest frames to the lab frame can be 
schematically described as: 

^CM^ ' ' Pl^ I I 

qi,uji > q^,u^ > q^,uj-^ 

(72, W2 > q2,^^2 ^92,^2, (3) 

where q{, (jj , and UJ2 are the lab frame four- vectors of 
Qi, Q2, Xi and X2, respectively. 

In practice, \Pt\ ~ pip^/\/l ~ pfp^'/2MG, where p^'^''^ 
is the magnitude of the vectorial sum of the transverse 
momenta of initial state radiation. As a result \Pt\ ^ 1 
for sufficiently large values of Mq. We will use the ap- 
proximation /3t — >■ for the remainder of the discussion. 

The relevant experimental observables are the three 
momenta of the massless final state particles, Qi and 
Q2 and the missing transverse momentum, denoted M , 
which is the vectorial sum of the transverse momenta of 
the particles xi and X2- Even with the aforementioned 
approximation Pt — >■ 0, it is not possible to fully solve 
for each of the kinematical unknowns using the available 
observables in this example. 

However there is an additional, well-motivated approx- 
imation that can be made to significantly simplify the 
problem. If the mass Mq is sufficiently large, relative to 
the hadron-hadron collider energy ^/s, the particles Gi 
and G2 will be mostly produced near the VI ^ 2Mg 
threshold, such that jcm ^ 1- 
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Assuming ^cm = 1 implies, of course, that /3cm = 0. 
In this approximation, the CM frame of this hard process 
is now also simultaneously the rest frame of the particles 
Gi and G2. The momenta of their decay products is 
given by Equation [2] with qi,ijJi = q^jUj^. Importantly 
1 91 1 — I (ft I — Ma/2 in this frame. 

Hence we can move from the laboratory frame to the 
CM frame by finding the longitudinal boost to a refer- 
ence frame where the magnitude of the momenta of the 
objects Qi and Q2 are equal. We denote this reference 
frame the rough approximation-iiame, or i?-frame, and 
the longitudinal boost moving from the lab frame to the 
i?-frame as The momenta of Qi and Q2 in the R- 
frame are denoted and respectively. With the 



constraint |gf | — we find that 
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(4) 



Furthermore, we define the i?-framc mass, M^, as 



A'^'qj^Mq and integrating over Xa we find that the differ- 
ential cross-section with respect to 7cm is given by: 
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[flMf2i 



d'ycM 
47, 



SXn 



S"fCM 



a ^ b] 



dXn 



(8) 



In Fig. [T] we show the probability distribution function 
for 7cAf for y/s = 14 TeV, where we have numerically 
integrated Eq. [Sjfor qq-\\ke {u and sea quark PDF's) and 
gg-\ike production. We use PDF parameterizations of the 
form xf,{x) = AiX^^{l-x)'J^{l + eiy/^ + -/ix) + A'iX^'^{l- 
x)^i with NNLO parameters determined from a global 
PDF fit at Q2 = 1 GeV^ [16]. Larger values of Mc result 
in lower values of 7c _a/, with all distributions peaking 
at approximately 7c j\/ ^ 1-1 and falling quickly with 
increasing 7cm- 



Mr = 2|qtl = m 



ill. - (fizY - W 



4o)^ 
(5) 

As 7cA/ 0, we find that q^ — >■ qi, fin — >■ (5i and 
Mr — > A/a. It should also be noted that the quantity Mr 
is invariant under longitudinal boosts so even if 7cm 7^ 1, 
Mr is independent of the true value of Pi . 



III. THE 7CM = 1 APPROXIMATION 

In order to understand how ^cM is distributed we con- 
sider the simple model with two scalar particles: $0 with 
zero mass and $1 with mass Mq. Using the notation of 
Sec.|lTj we consider contact interaction pair production of 
<I>i through a $0$^ vertex. The sub-process cross-section 
is proportional to 



ays) cx A ; — - — oc A 



7cmM 



(6) 



where A is the dimensionless <I>g$^ coupling, which we 
set to 1. From Eq. |6] we observe that 7cm = 1 is kine- 
matically forbidden, and that the cross-section for the 
sub-process will decrease asymptotically as 1/7ca//- 

Additional suppression of large values of 7c m is caused 
by the parton distribution functions (PDFs) in hadron- 
hadron collisions. Assuming the two initial state 4>o par- 
ticles are partons from colliding protons with momentum 
fractions Xa and Xh respectively, and PDFs fi(x) and 
f2{x) we can write the total cross section as: 



da 



dXadXb 



« [fliXa)f2{xb) + a O b]a{s = SXaXb), (7) 



where s is the proton-proton CM energy. Changing vari- 
ables from Xb to 7cM through the relation sXaXb = 
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FIG. 1: Distribution of 7c a/ for qq-like and gg-\ike production 
at ys = 14 TeV for different values of Mq. 



The exact dependence of the sub-process cross-section 
on 7cM will vary depending on the nature of the in- 
teracting final and initial state particles in the 2 — !• 2 
process, but the resulting distribution of 7c a/ should be 
qualitatively similar to the result show in Fig. [T] 7c Af 
exactly equal to 1 is kinematically forbidden, but values 
of 7cAf near 1 are preferred to larger values due to the 
falling sub-process and total cross-sections with increas- 
ing \/l = ^SXaXb oc 7CAf- 

To derive the expression for Mr in Eq. [sj we used the 
approximation 7c a/ precisely equal to 1, and found that 
Mr — >• Ma as 7c a/ — 1. In order to understand the be- 
havior of Mr when 7c a/ 7^ 1, we return to the example 
introduced in Sec. |llj Using the same notation, we again 
consider the pair production of massive particles Gi and 
G2, and continue to use the approximation (3t 0, this 
time with jcm not equal to 1. We recall that the vari- 
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FIG. 2; Distribution of Mr, in units oi '^cmMa, for different 
values of 'ycM- Distributions are normalized such that their 
maximum value is equal to one. 



able Mr is invariant under longitudinal boosts, so its 
value is independent of the true value of /3; . This implies 
that for fixed 7cm, Mq and M~^, there are six remain- 
ing degrees of freedom, corresponding to the unit vectors 
PcMi Ui and 1*2- We numerically integrate over these 
angles, making the assumption that each of these unit 
vectors are independent of each other with flat probabil- 
ity distribution functions on the unit sphere, as if all the 
interacting particles are scalars. The resulting distribu- 
tions of Mfl, for different fixed values of 7cm, are shown 
in Fig. [2j We observe that the peak value of Mn scales 
as ICM^Ai with the width of the M/j distribution in- 
creasing with ^CM- Hence, in practice, the distribution 
of Mr will peak near Ma, even when 7c a/ 7^ 1- 



IV. SIGNAL AND BACKGROUND 
DISCRIMINATION: THE RAZOR 



In Sec. Ill we found that the distribution of the observ- 
able Mji contains information about the mass difference 
Ma, introduced in the example described in Sec.|lTj This 
is a useful property for characterizing this process. We 
further explore whether and when this variable is also 
useful for selecting events in the presence of background, 
establishing the discovery of this process. This depends 
on the experimental final state considered and the signal 
and relevant backgrounds. 

We consider massive di-squark production in a generic 
supersymmetric extension of the Standard Model, where 
each squark decays directly to the lightest supersym- 
metric particle (LSP) and a quark, with the LSPs es- 
caping detection. The relevant final state contains two 
or more jets and missing transverse energy. A particu- 
larly challenging background to this final state is QCD 



multijet production, where non-zero missing transverse 
energy can result from instrumental backgrounds, jet 
mis-measurements, finite detector acceptance and non- 
gaussian tails in the detector response, in addition to the 
production of neutrinos within jets. 

Using the notation of Sec. [IT] to describe our signal 
process, the particles Gi and G2 are the squarks with 
mass Mq = Mq, Qi and Q2 are the quark jets, and xi 
and X2 are the LSPs with mass M^. The final state 
observables are the two jet four- vectors, q{ and (73, and 
the missing transverse momentum, M. 

For the background processes, we consider QCD dijet 
production. In the dijet rest frame, we express the two 
jets' four- vectors as: 



^2 = —{1,-1)}, 



(9) 



where Vs is the dijet invariant mass and u is a unit vec- 
tor. If we assume that the Lorentz transformation from 
the dijet rest frame to the laboratory frame is simply a 
longitudinal boost, /3/, we find that for this type of event 
Mr — Vs. Therefore, Mr will be distributed as VJ for 
this background process, falling steeply, while the signal 
distribution will peak near M^ . The question of whether 
or not we can identify signal events in the presence of this 
background becomes a question of whether the effective 
dijet cross-section is sufficiently small for in the range 
of the signal peak around Ma , given a set of event selec- 
tion requirements. 

Traditionally, the requirements that have been used to 
improve the signal to background ratio in such processes 
is large magnitude of the missing transverse energy, M, 
as well as large transverse momenta of the two jets. These 
are highly correlated with the observable Mr. For the 
background process we are considering, the magnitude of 
the missing transverse energy (from mis-measurements of 
the dijets' transverse momenta or neutrinos produced in 
the jets) is highly correlated with the magnitude of the 
dijets' transverse momentum which, in turn scales with 
\/§. These observables are also highly correlated for the 
type of signal events we are considering since the LSPs' 
transverse momenta will also scale with Ma, and as a re- 
sult so will the missing transverse energy. A requirement 
on each of these variables is a requirement on the scale 
of the signal and background events. As a result, for a 
given integrated luminosity, if the background yield in 
the region \/~§ ^ AI^ is prohibitively large relative to the 
signal yield, it is unlikely that additional hard require- 
ments on the magnitude of the missing transverse energy 
or the jet transverse momenta will assist the discovery. 

Examining the expression for Mr in Eq. [5] we see that 
there is additional kinematical information not yet used. 
For example, Mr is independent of the azimuthal angle, 
A0, between the two final state jets. For the QCD dijet 
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background, the jets should be largely back-to-back in 
the transverse plane, with A(/) peaking at tt. On the other 
hand, the two jets in the SUSY signal events result from 
the decay of two separate squarks, implying that their 
direction in the transverse plane is largely independent 
of each other, apart from spin-correlations and effects 
resulting from PcM 7^ 0. Hence, the distribution of 
for signal events will be significantly flatter than for the 
background. Rather than simply cutting on the variable 
A</), we incorporate this information into a new variable 
denoted . 

In this particular final state, we assume that there 
are two escaping weakly interacting particles, whose 
four-momenta in the laboratory frame we denote i/J 



and 



^2; 



with each particle "paired" with an ob- 



served jet with four- momenta q\ and ql^: respectively. 
Fr om these four-vectors we de fine the variable M2G = 
^/{l/2)[{v{ + q[y + (i//, + which is equal to for 

signal events. The only constraint we have on the four- 
vectors i^l is that the vectorial sum of their transverse 
momenta should be equal to the observed missing trans- 
verse energy, M. Setting (i^')^ = and minimizing M2G 
over and i^2z yields: 

niinAf2G = I^itII'^itI ^ 9it ' ^ir 

(10) 

Motivated by the background we are considering, we 
assign half of the measured missing transverse momenta 
to each escaping particle so that i7[rp — v\rp — M /2 and 
define as: 



92T ' '^2T 



M? = y ^(klrl + klrl) - \m ■ (91t + ^It) ■ (H) 

The variable also contains information about the 
scale of the process we are studying. If we assume that 
7cm = 1 then the Mi^ distribution has an endpoint 
at Ma for signal events. We note that is an ad- 
ditional measurement of the scale of the process that 
uses information independent of the M^. Therefore, 
rather than cutting on we form the dimension-less 
'i?-frame razor', i?, as the ratio of M|? and M/j, such 
that R = M^/Mr. For the signal process, the distribu- 
tion of R peaks near 0.5, since this is the ratio of two 
measurements of the same scale, A/a, with an additional 
geometrical factor due to the fact that Af|? contains only 
transverse information. For the QCD dijet background, 
if A? = 0, then i? is 0, for any value of v^. 

As was discussed previously, there are several mech- 
anisms for the measurement of M to be non-zero in 
QCD dijet events. For example, one or both jets in 
the final state could be mis-measured due to calorime- 
ter non-compensation, uninstrumented regions of the de- 
tector or weakly interacting particles produced within 
the jets, causing an imbalance in the event and result- 
ing in non-zero missing transverse momentum. To evalu- 
ate how these possibilities affect the measured values for 



Mr and in background events, we return to Eq. [o] 
which describes the kinematics of the dijet system in it's 
CM frame. We now assume that the measured jet mo- 
menta, q\, are scaled relative to their true values, so that 
q\ — fiki. Here, we are assuming that the direction of 
the two jets is not changed, but rather that only a frac- 
tion fi of the jets' momentum is observed, where /j > 0. 
Additionally, without loss of generality we adopt the con- 
vention /l > /2. 

With these mis-measurements, we find that Mr takes 
a value: 



Mr, 



(12) 



independent of the longitudinal boost, that takes the 
jets from their CM frame to the laboratory frame. The 
missing transverse energy can now be non-zero, with 
M = (/2 — fi)kiT and M^ can be expressed as: 



= \/(/l-/2)/l 



s(l — {v ■ z)^) 



(13) 



From Eq. 12 we see that these mismeasurements de- 
crease the value of Mr, assuming that /i ^ 1. There- 
fore the distribution of Mr for the background will not 
have events promoted to the tail of the distribution due 
to these types of mis-measurements; instead, these mis- 
measurements will suppress the background Mr distri- 
bution. Furthermore, if we require that R > C, where 
C is some cut value, this implies that CMr < M^j}. To 
understand the effect of this cut, we change variables 
{v ■ zf = cos(0i)2 and /i = /2Cos(6'2)^ With these 
substitutions, we re-express the inequality CMr < Af^ 
as: 



16C^ cos(6'i)2cos(6l2)'^+ sin(6'i)Sin(6'2)^ < 
4sin(6'i)2 cos(0i)2 sin(6l2)^ cos(6i2)^. 



(14) 



This inequality implies that if C > 1/2, no background 
events of this type will satisfy the requirement on R. if 
C ~ 0.4, some events can pass, but M^ will reach its al- 
lowed maximum, for fixed -s/I, at M^ ~ with the 
razor inequality implying that Mr < M^/C < Vs/2. 
Hence for this type of background event to result in 
Mr ~ Ma, it must have > 2Ma- Therefore, we 
observe that adding a requirement on R to our event 
selection will remove most QCD dijet events with mis- 
measurements of the type described above. 

Another possibility resulting in non-zero missing trans- 
verse momentum in these background events is that there 
are additional particles, whose vectorial sum of trans- 
verse momentum is non-zero, that escape detection. For 
example, jets resulting from initial state radiation could 
remain unseen due to limited detector acceptance, caus- 
ing a transverse imbalance in the visible momentum in 
the event. In order to understand the effect of this type 
of background on Mr and M^, we consider the follow- 
ing example. We denote the vectorial sum of the trans- 
verse momentum of particles escaping detection as Pt- 
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Returning again to the QCD dijet example described 
by Eq. ^ a nonzero value of Pt will result in two sig- 
nificant changes to the final state particle kinematics. 
Firstly, the missing transverse energy will be non-zero, 
with M = Pt- Secondly, this missing momentum will 
result in the dijet system undergoing an additional trans- 
verse boost when moving from the dijet rest frame to the 
laboratory frame (any additional contribution to the lon- 
gitudinal momentum imbalance in the event is absorbed 
into the longitudinal boost, /3/, which moves the dijets 
from their CM frame to the laboratory frame). Specifi- 
cally, the dijets are moved to a frame traveling at a ve- 
locity (3 = M/{jV§), where 7 = (1 - |/3p)i/2 and VS is 
the dijet invariant mass. In this case, Mr is given by 



-1/2 



(15) 



while M^ can be expressed as 



with equal mass, that both decay directly to a mass-less 
visible particle and a massive invisible particle. Despite 
the fact that Mr is motivated by this particular type 
of example, we find that it is useful in a more general 
context. 

We return to the example described in Sec. |ll] but we 
now allow for the two massive particles, Gi and G2, to 
have different masses. Alternatively, we can assume that 
the weakly interacting particles resulting from the decays 
of Gi and G2 have different masses. Using the notation 
of Sec. HD we will assume that each of the two decay 

, so that 



-MZ 



chains has a different value for Af^ 

Ml^MX{l + S)=MAil + S). 

Assuming 7c a/ = 1, we numerically integrate over the 
angular degrees of freedom contained in the variables 
Ml, U2 to derive the distribution for Mr^ for different val- 
ues of 6, which is shown in Fig. [3] We find that Mr peaks 
precisely at the geometric mean of and M^. 



M^ : 



/7/3s(V(l-({)-z)2+7/3) 



(16) 



We observe that that for fixed vl, after applying a re- 
quirement on R, remaining background events will have 
Mr with an upper bound that goes as if the jets 

have a large transverse component in their rest frame, 
otherwise as 7/3\/i. Recalling that 7/? — \Pt\/\/~^, we ob- 
serve that the asymptotic behavior of these upper bounds 
can be re-expressed as \Pt\ and {\Pt\V^Y^'^ ^ respectively. 
Hence, we see that in order for these types of background 
events to populate the Mr distribution in the neighbor- 
hood of some value of M\ , the magnitude of the vectorial 
sum of the transverse momentum of any missing particles 
needs to be on the order of Ma. 

In the case of the jets plus missing transverse momen- 
tum final state, this example is not only relevant for the 
QCD multi-jet background, but also for the so-called ir- 
reducible background Z {1/1/)+ dijets. Here, \Pt\ ^ , 
and hence has an intrinsic scale on the order of Mz- The 
distribution of Mr, still falls off exponentially for this 
background when Mz < Mr. 

We observe that Mr is potentially a powerful variable 
for distinguishing SUSY dijet plus missing transverse mo- 
mentum events from the relevant backgrounds when used 
in conjunction with requirements on the i?-frame razor, 
R. In particular, we note that this variable is robust 
against effects related to jet mis-measurements and lim- 
ited detector phase-space acceptance which often result 
in spurious missing transverse momentum. 



V. GENERALIZING THE APPLICATION OF 

Mr 

In Sec. |Tl]we introduced and derived the variable Mr 
in the context of pair production of two massive particles. 
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FIG. 3: Distribution of Mr, in units of MaVI + S, for dif- 
ferent values of 5. Distributions are normalized such that the 
maximum value is equal to 1. 



To define Mr in cases with more than two visible par- 
ticles in the final state, we generalize the two-object to 
a multi-object final state by forming two pseudo-objects. 
Ha and Hi,, with four- vectors ha and hb, respectively. 
Each pseudo-object's four-momenta is simply the sum of 
the four- vectors associated with it. Each final state ob- 
ject can only be assigned to one pseudo-object, and each 
pseudo-object must have at least one object assigned to 
it. These associations are determined by minimizing the 
quantity (/la)^ + (^b)^- Once the two pseudo-objects are 
defined, their masses are set to zero with the direction 
and magnitude of their momenta unchanged. With this 
prescription, AIr is defined as in Eq. [5j with q{ = ha and 

For example, we assume that the particles Gi and G2 
have the same mass {Mq), as do xi a-^d X2 (M^), except 
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now one or both of the particles Gi undergoes a two-body 
decay to a visible particle, Qi+2, and another particle, Si, 
with mass Ms = Mq{1 — 5). The particle Si then de- 
cays to another visible particle, Qi, and Xi- Numerically 
integrating over all the decay angles in this scenario (us- 
ing angular probability distribution functions flat on the 
unit sphere) with 7cm = Ij and requiring R > 0.4, we 
derive the distributions for M/j, assuming one or both of 
the particles Gi decays to an intermediate Si, shown in 
Fig. |4] We find that, in both of these cases, the result- 

M^ —M^ 

ing Mfl, distribution peaks at Ma = ^ , regardless 



the scale 



Mi -Mi 



of the value of S (for the values considered here) and ir- 
respective of whether all of the visible decay products 
resulting from a particular Gi are assigned to the same 
pseudo-object. 
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FIG. 4: Distribution of Afjj when one (left) or both (right) 
of the particles d decays to an intermediate particle Si with 
mass Ms = Mg{1 — 5), for different values of S. Distributions 
are normalized such that the maximum value is equal to 1 



A. Example: inclusive search for SUSY 

A potential use-case for an analysis incorporating the 
variables Mfj and R is an inclusive search for R-parity 
conserving SUSY. Generally, these scenarios can be char- 
acterized by the production of pairs of massive, strongly- 
interacting super-partners decaying directly or through 
cascade decays to SM particles, and LSPs that escape 
the detector unseen. 

The canonical example described in Sec. [IT] is the sim- 
plest case of this type of SUSY particle pair production, 
with the particles Gi representing squarks that each de- 
cay directly to a quark and neutralino. The example 
introduced earlier in Sec. |V] of two heavy particles un- 
dergoing two-body decays to intermediate heavy parti- 
cles, which subsequently decay to a visible and weakly- 
interacting particle, describes the pair-production of 
gluinos which decay to a quark and squark, subsequently 
decaying as in the canonical example. This cascade decay 
is often quite prominent when Mg > Mq. 

Hence, for R-parity conserving SUSY scenarios, di- 
squark production will result in an Mn peak around 



with potentially several dif- 
ferent peaks corresponding to the different squark gener- 
ations that can only be resolved if the masses are suf- 

> Mfj, gg and gq produc- 



ficiently different. If Mg > Mq, 



tion will result in M^ peaks at 



M5 



and 



Mi-Mi 
Ms 



The result is Af^ spec- 



Ml^ = 

troscopy, with with different peaks corresponding to the 
respective mass differences between the massive SUSY 
particles produced first in the hard scattering process 
and the LSP, and the geometric means of these mass 
differences when two SUSY particles are produced with 
different masses. 

If one approaches the jets -I- missing transverse en- 
ergy -|- X final state in an inclusive way, the dominant 
backgrounds will be QCD multi-jets, tt-|-jets and y-|-jets, 
where ^ is a IF' or Z vector boson decaying to leptons 
and/or neutrinos. The QCD contribution to this back- 
ground will be largely marginalized by the use of a cut 



on the razor, R, as described in Sec. IV The A/^ distri- 
bution of QCD events passing this requirement will fall 
roughly exponentially in Af/j, with a slope largely deter- 
mined by the value of the cut on R, but not exceeding 
the slope of the for these QCD processes. The re- 
maining backgrounds identified above must have a large 
transverse momentum imbalance in the event in order to 
pass the requirement on R, implying that their contri- 
bution to the distribution in Mji will be comprised of 
events with final state neutrinos or leptons that are not 
explicitly identified as such. As was discussed in Sec. |IV| 
with the example of Z(j/i/)-f jets, the Mn distributions 
for these processes will also fall roughly exponentially, 
with slope determined predominantly by each processes' 
respective scale, Mz, Mw and Mj. 

As a result, if the various values of Ma that charac- 
terizes a particular SUSY scenario are sufficiently higher 
than Mz, Mw and Aft, the signal events will appear as a 
wide peak(s) on top of falling exponential backgrounds, 
potentially a striking signature depending on the relative 
production rates of the SUSY and background processes. 



B. Example: H W{iu)W{lv) 

We have, so far, only discussed the canonical scenario 
described in Sec.|lT]in the context of cases where 7cm is 
near one, corresponding to the pair-production of mas- 
sive particles near threshold. The same example ap- 
plies to the process of a Higgs boson decaying to two 
IF' bosons, which subsequently decay leptonically. Now, 
ICM = Mh/2Mw 

As we saw in Sec. |III[ the peak value of Mu will scale 
as jcmMa which, in this case, implies that Mr will peak 
roughly at Mh/2. Fig. [5] shows the distribution of Mji 
for different values Mh, with 

In practice we know that 7^ 0, and, in fact, it scales 
with Mh- In Fig. |6]we illustrate the effects of non-zero 
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2Mr / Mjj 



FIG. 5: Distribution of M_r for the di-lepton final state in 
H — >■ WW as a function of M_rf, where we have made the 
approximations = and neglected spin-correlations. 



pf[ . We show in Fig. [g] (left) the distribution of for 
different values of Mh, with pi[ = Mh/5. In Fig. [g] 
(right) we illustrate the situation where Mh = 170 GeV 
for different values of . 
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FIG. 6: Distribution of Mr when: (left) = Mr/5, (right) 
AIh ~ 170 GeV. Spin-correlations are neglected in both cases. 



We observe that 2Mj^ peaks around the true value of 
Mh- This peak is naturally quite wide, such that the lep- 
ton reconstruction resolution will have a negligible effect 
on its shape. An analysis using AIn will also benefit from 
a selection including a requirement on R, which should 
marginalize backgrounds where one or both leptons re- 
sults from weakly decaying hadrons. Other backgrounds, 
such as Z{££) and tt production, as was the case in the 
jets -|- missing transverse energy final state described in 
Sec. |IV[ will fall roughly exponentially in AIjj once M/j 
exceeds their respective mass scales. 



VI. Mr, 

In Sec. [nj we described how to move from the labora- 
tory frame to the i?-frame via a longitudinal boost, /3_r. 
With the assumptions described in that section, partic- 
ularly that 7cM = 1, we found that the variable 
is equal to Ma- In subsequent sections we discussed 
the properties of the variable when some of these 
assumptions no longer apply and for deviations from 
the simple scenario described in Sec. [llj We have shown 
that the useful properties of the variable Mji are robust 
against the variations we have considered, but there is an 
important caveat. 

Namely, when 7c Af deviates from one, there are situ- 
ations when > 1, such that no longer describes 
a physical boost and the i?-frame is ill-defined, as is the 
variable Mr. In this section we will describe a set of ref- 
erence frames, which we will denote the i?*-frames, which 
are always well-defined, and variables R* and M/j. . 

Firstly, we identify the kinematical characteristics that 
are associated with events with \f3fj\ > 1. Using the no- 
tation of Sec. [nJ we again consider the pair-production 
of massive particles resulting in two visible particles, Qi 
and Q2, along with missing transverse energy M. Set- 
ting 7c Af — 1-1, we scan over values for the unit vectors 
Ml, M2 and PcM, noting for which values and with what 
frequency we find > 1. In Fig. [T] we show the corre- 
lation between the normalized ^-components of momenta 
of Qi and Q2 in the rest frames of their respective par- 
ents Gi for events where the i?-frame is ill-defined. We 

find, as perhaps one could infer from the expression of 

I _ I 

Pr = — that these longitudinal momentum com- 
ponents tend to be equal in both direction and magni- 
tude. In fact, as ^cm tends toward one, the distribution 
shown in Fig. [7] tends toward a discrete line along the 
ill ■ z = U2 ■ z diagonal. 

In Fig. [8] we show the correlation between the difference 
in azimuthal angles between the momenta of Qi and Q2 
and between the momenta of Qi and Pcm- We find that 
events with > 1 tend to have ui and U2 pointing 
in the same direction in the transverse plane, with Pcm 
pointing in either the same or opposite direction. 

These observations indicate that the cases where the 
i?-frame is ill-defined result from the neglecting of the 
transverse component of j3cM in the approximations 
made in the derivation of Mn^,. Here, we describe the 
derivation of variables that take into account this trans- 
verse component. 

We consider a longitudinal boost to the four-vectors 
q{ and associated with a velocity /3l*- Subsequently, 
we consider a transverse boost, with velocity , that 
is applied in opposite directions to the four-vectors as- 
sociated with Qi and Q2, taking them to different ref- 
erence frames. We denote these two reference frames as 
i?*-frames, with the additional requirement that, in each 
of their respective i?*-frames, the momenta of Qi and 
Q2 must be equal in magnitude. The series of Lorentz 
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FIG. 7; Correlation between Ui ■ z and 112 ■ z for events with 
7CM = 1.1 and \I3r\ > 1. Distribution is normalized to unit 
volume. 




A(|)(U^,U2) 



FIG. 8: Correlation between A4){ui,U2) and A(j){ui, /3cai) 
for events with 7c a/ ~ 1.1 and > 1. Distribution is 

normalized to unit volume. 



which can be used to solve for |/3^ | = in terms of 
Pf and 



/3#' ■ (g*iT + 9*2t) 



(19) 



Just as we did in the i?-frame, we will define the R*- 
frame mass, M^* , as two times the magnitude of the 
momentum of Qi in it's respective i?*-frame. Af^j. can 
be expressed as: 

Mfl. ^2|gfV2kTl = 



(91 



IT 



920 



(20) 



In order to calculate Mjj* , we must specify values 
of 13 and . Motivated by the event configura- 
tions which lead to > 1, in particular those de- 
scribed in Fig. pi we choose a value of (an angle 
in the azimuthal plane) which maximizes the quantity 

l/^T ' ('?iT + 92t)I- With this choice, (3^ can be ex- 
pressed as: 



3^' 



9lT 92T 



(21) 



The value of /3l* is determined by requiring that the 



condition 



is satisfied. This choice results in 



/3i 



4iz 



(22) 



With each unknown quantity now specified, Affl. can be 
expressed, event-by-event, as 



|9it + '?2tP 

(23) 

Another interesting quantity is 7ij* =(1 — P)~^^^, 
which can be expressed in terms of lab frame observables 

as 



boosts taking Qi and Q2 from the laboratory frame to 
their respective i?*-frames can be summarized as: 



q[^q1^ qf 
q'2^q'2^ qf 



(17) 



The constraint that Qi and Q2 have the same energy 
in the respective R* frames can be re-expressed as a con- 
straint equation on the variables and : 

lL'iq[o-qw)-lL'l3L'iqL-qL) = A? ' ' (^It+^It) (18) 



IR' = 



iq[o + q2o)^ - (qL + qL)^ 



\ iq{o + qioy-i<iL + qLy 



{24) 



As is the case for Af^,, Mfj* is invariant under longitudi- 
nal boosts, as is 7^* . 

Similarly as for the i?-frame, we define the i?*-frame 
razor, R* , as the ratio of M^' and Af^. , with AI^ given 
by Eq.[Tl] 

To understand how the distribution of Mr* changes 
with 7cM, we numerically integrate over all the decay 
angles, assuming their distributions are flat on the unit 
sphere. The resulting Af andjRtMfj* distributions are 
shown in Fig. [9j We find that the peak value of the Mr* 
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FIG. 9: Distribution of Mr* (left) and ^r-Mr* (right) for 
different values of 7c a/- Distributions are normalized such 
that their maximum value is equal to one. 



distribution is at approximately Ma, regardless oi^cM^ 
while ^R*Mii* peaks at jcmM^. 

Comparing Fig. [9] and Fig. [2] we see that the peak 
position of the jr* Mr* distribution scales like the peak 
of the Alfj distribution. Mr is a variable most useful for 
treating the case 7cm = 1 which, in practice, is kinemat- 
ically forbidden. The quantity ^r^Mr', reproduces the 
same peaking behavior, without ill-defined configurations 
and better resolution on the quantity JcmM^. 

In fact, the variables Mr-, jr-Mr* and Mr share 
many properties. We consider two of the examples from 
Sec. |V) now in the context of Mr* and jr> . The first 
scenario is of two massive particles, Gi and G2, with 
different masses decaying each to a visible particle and 
potentially massive weakly interacting particle, such that 
Ml = Mi(l + S) = Ma(1 + S). Assuming 7cm = 1, 
and numerically integrating over the angular degrees of 
freedom assuming scalar decays, we calculate Mr* as a 
function of S, with the resulting distributions shown in 
Fig. 10 (left). We observe that Mr* , like Mr, has a peak 
whose position scales with \/l + S. 
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-5=0.3 
■■ 5 = 0.5 
-■5= 1.0 




M„., / M\' 1+S 



0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



FIG. 10: (left) Distribution of Mr* , in units of MaVI + 5, for 
different values of S. (right) Distribution of 7fl« Mr* when one 
of the particles d decays to an intermediate particle Si with 
mass Ms = Mg(1 — 5), for different values of 5. Distributions 
are normalized such that their maximum value is equal to one. 



The second example from Sec. [V| involves two particles 
Gi and G2, with the same mass. Gi undergoes a two- 



body decay to a visible particle, Q3, and another particle. 
Si, with mass Ms = Mc{l — S). The particles G2 and 
Si then each decay to a weakly interacting particle and a 
visible particle, where the mass of the weakly interacting 
particles is Af^ . The numerically integrated ^r* Mr* dis- 
tributions, for different values of 5, are shown in Fig. [lO| 
(right). We observe that, like Mr, the quantity ^r*Mr* 

— M^ 

peaks at Ma = ^ , regardless of the value of 6. 

The interplay between R* and Mr* {'Jr* Mr* ) is qual- 
itatively the same as between R and Mr for di-jet back- 
grounds, as described in Sec. |IV[ A requirement on R* 
suppresses contributions from mis-measured di-jet events 
significantly, and the quantity 7^* Mr* peaks at The 
asymptotic behavior of Mr* is similar for backgrounds 
with an escaping (weakly interacting or outside accep- 
tance), high transverse momentum, particle or system 
of particles. We note that Mr* {jr* Mr* ) can be used 
similarly to Mr in the context of an inclusive SUSY 
search, as described in Sec. |V A| or in the search for 
H WW Mv, as described in Sec. [VB 



In the latter case, we have an additional piece of in- 
formation when using Mr* {"/r* Mr* ) rather than Mr 
in 7^'. We find that for H -)■ WW lv(.v decays, 
Mr* wiU peak at Mw while -i^' Mr* will peak at Mh/1. 
These two observations are demonstrated in Fig. |11| 
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FIG. 11: (left) Distributions oi Mr* and (right) ')r*Mr* for 
H — > WW — >■ Ivlu events.. We make the approximation that 
Pt — and spin correlations are neglected. 



A. Outlook 

We introduce a set of variables Mr and Mr* designed 
to study the characteristics of processes involving the 
pair-production of massive particles that each decay 
directly or through a cascade of decays to SM particles 
and weakly-interacting particles escaping detection. We 
also describe the dimension-less variables R and R* 
which, used in conjunction with Mr and Mr* , provide 
the means to select these processes in the presence of 
large backgrounds, for a variety of final states. We have 
described how these variables can be used to discover 
and characterize R-parity conserving SUSY scenarios 
and SM Higgs boson decays to leptonically decaying 



10 



VF- bosons. 



hosting a talk on the subject. 
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